Quasienergy description of the driven Jaynes-Cummings model 
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We analyze the driven resonantly coupled Jaynes-Cummings model in terms of a quasienergy 
approach by switching to a frame rotating with the external modulation frequency and by using the 
dressed atom picture. A quasienergy surface in phase space emerges whose level spacing is governed 
by a rescaled effective Planck constant. Moreover, the well-known multiphoton transitions can be 
reinterpreted as resonant tunneling transitions from the local maximum of the quasienergy surface. 
Most importantly, the driving defines a quasienergy well which is nonperturbative in nature. The 
quantum mechanical quasienergy state localized at its bottom is squeezed. In the Purcell limited 
regime, the potential well is metastable and the effective local temperature close to its minimum 
is uniquely determined by the squeezing factor. The activation occurs in this case via dressed spin 
flip transitions rather than via quantum activation as in other driven nonlinear quantum systems 
such as the quantum Duffing oscillator. The local maximum is in general stable. However, in 
presence of resonant coherent or dissipative tunneling transitions the system can escape from it 
and a stationary state arises as a statistical mixture of quasienergy states being localized in the two 
basins of attraction. This gives rise to a resonant or an antiresonant nonlinear response of the cavity 
at multiphoton transitions. The model finds direct application in recent experiments with a driven 
superconducting circuit QED setup. 

PACS numbers: 78.47.-p, 74.50.+r, 42.50.Pq, 42.50.Hz 



I. INTRODUCTION 

Damped nonlinear classical oscillators can display 
rather nontrivial features when they are modulated by an 
external time-dependent driving forced. The simplest ex- 
ample is the classical Duffing oscillator, for which a quar- 
tic nonlinearity extends the harmonic potential such that 
the static potential still remains monostable. Adding a 
periodic modulation generates several dynamical stable 
states and the role of environmental fluctuating forces be- 
comes particularly intriguing at the bifurcation points^. 
Close to the fundamental resonance, the classical Duffing 
oscillator displays two stable states of large and small os- 
cillation amplitudes. Environmental fluctuations can in- 
duce transitions between the two states and the scaling 
property of the probability of an activated escape from 
such a metastable state near the bifurcation point can be 
determined^i^. 

For the corresponding quantum Duffing oscillator, two 
inherently quantum mechanical effects can contribute to 
the escape from the metastable forced state over the ef- 
fective quasienergy barrier. One mechanism is dynam- 
ical tunneling^ in the quasienergy surface, which leads 
at low temperatures to a sharp increase of the transition 
probabilities near the classical bifurcation point (deter- 
mined in the quasiclassical approximation). In addition, 
quantum activation^ occurs which even at zero temper- 
ature leads to an effective diffusion over the quasienergy 
barrier induced by the environmental quantum fluctua- 
tions. At a flrst instance, this effect might appear some- 
what counterintuitive, but one has to bear in mind that 
energy absorption in the environment can increase the 
system's quasienergy, this quantity being defined in a 



rotating frame. It has been shown that the escape is al- 
ways over the barrier (of activation type) , provided that 
the broadening of the quasienergy levels, induced by the 
environment, exceeds the corresponding coherent tunnel- 
ing rate^. Other driven nonlinear quantum systems, such 
as the parametrically driven Duffing oscillator— or mod- 
ulated large-spin systems^ show quantum activated be- 
havior as well. In this regime, there is a separation of 
time scales. After the fast intrawcll relaxation processes 
have occurcd, the system occupies a quasiclassical state, 
which is metastable. From this, it can escape by rare 
interwell transitions that eventually would lead to the 
other metastable state. 

In the opposite limit, when tunneling transitions no- 
ticeably contribute, the separation of time scales is not 
well defined and a mixture of all quasienergy states forms 
the stationary state. This is the case when a coher- 
ent resonant excitation induces a population of a mul- 
tiphoton quasienergy state, as shown for the Duffing 
oscillator—"—. Since the quasienergy states oscillate with 
different phases with respect to the external modulation, 
a resonant or antiresonant response of the oscillator to 
the modulation may occur at a multiphoton transition. 
Which type arises, is determined by the dominant sta- 
tionary population of the involved quasienergy state and 
thus depends also on the parameters of the environment. 
These lineshape properties around a multiphoton (anti- 
)resonance are connected with a resonantly enhanced es- 
cape in form of resonant dynamical tunneling^ii^, which 
shows up as resonant tunneling peak in the switching 
rate. Enhanced peaks are associated with resonant line 
shapes, while reduced peaks go with an antiresonance. 

Recently, we have showni^ that a dynamical bistabil- 



ity also occurs in the setup of a driven linear resonator 
coupled to a quantum two-level system (driven circuit 
quantum electrodynamics (QED) set-up). The system 
is conveniently modeled by the Jaynes-Cummings (JC) 
modeli^ which is extended by a driving term. The static 
Jaynes-Cummings model was originally studied to de- 
scribe the interaction of a two-level atom and a single 
quantized electromagnetic field mode and has an inher- 
ent nonlinearity since the splitting of the vacuum Rabi 
resonance depends on the number N of photons in the 
resonator as vW+l. Already its undriven dynamics has 
many interesting facets, including Rabi oscillations, col- 
lapses and revivals of quantum states, quantum squeez- 
ing, quantum entanglement, Schrodinger cat and Fock 
states, and photon antibunchingi^. It also finds ap- 
plications beyond quantum optical set-ups, namely in 
nanocircuit architectures, such as Cooper pair boxes^^, 
superconducting flux qubita^^, Joscphson junctionsi^, 
and semiconductor quantum dots^^. In particular, the 
latter setups allow to explore the regime of strong qubit- 
resonator coupling as the resonator is typically formed 
by a nanoscalc on-chip transmission line. In addition, 
strong driving allows to access the regime of nonlinear 
response. 

These experiments based on quantum state engineer- 
ing with superconducting circuits were accompanied by 
progress in theory, based on the adoption of the undriven 
JC modeli^Ti^i to the particular experimental situations. 
In addition, the model has been extended by an addi- 
tional time-dependent modulation ter m^^i^'^ , mimicking 
the effect of rf fields applied either to the qubit or to 
the oscillator part. The more general case of N two- level 
atoms strongly coupled to a driven optical cavity has 
been considered in Ref. [2J for the case of weak driving 
fields and non-classical features in the photon correlation 
function have been identified. 

The driven Jaynes-Cummings model has recently 
found an experimental realization in form of a super- 
conducting transmon qubit device^^. The transmitted 
heterodyne signal has been used as a measure of the 
amplitude of the stationary oscillations of the oscillator 
position. It has allowed to study for the first time the 
supersplitting of each vacuum Rabi peak which reflects 
transitions between the groundstate and the first/the sec- 
ond excited state of the undriven Jaynes-Cummings spec- 
trum. Their energy difference is 2hg, where g is the in- 
teraction strength of the qubit and the harmonic mode. 
In addition, for stronger driving, the excitation of dis- 
crete multiphoton transitions up to a photon number of 
N = 5 has been observed. The measurements have been 
analyzed by accurate numerical simulations^. An effec- 
tive two-level approach which involves the vacuum and 
the one-photon dressed state has been used to describe 
the vacuum Rabi splittin g^^'^^ . Moreover, the character- 
istic vN spacing of the involved energy levels has been 
demonstrated. At resonance, the system is coherently 
excited to a iV— photon state. The relaxation occurs via 
subsequent dissipative transitions, generating eventually 



a steady state which involves a mixture of many quantum 
states and which renders a two-level description inappro- 
priate. 

In this paper, we complement our brief account of Ref. 
Il2l by a comprehensive theoretical analysis of the driven 
Jaynes-Cummings model in the deep quantum regime of 
few photons in the resonator. Motivated by the failure 
of perturbation theory for arbitrary small intensity of 
the driving, we carry on an alternative approach. The 
transformation to the frame rotating with the frequency 
of the external driving yields to a description in terms of 
quasienergy levels and states. Discrete multiphoton tran- 
sitions arise at avoided quasisenergy level crossings. Fur- 
thermore, we obtain the Hamiltonian in the dressed state 
basis in form of two quasienergy surfaces, one of which 
shows bistability. This description allows to interpret the 
multiphoton transitions as tunneling transitions in the 
bistable quasienergy surface. On the other hand, the re- 
gion of the spectrum for which perturbation theory fails 
can be studied by means of a harmonic expansion of the 
quasienergy surface around its global minimum. As al- 
ready mentioned in Ref. |l2i the lowest quasienergy state 
is a quantum squeezed state with the squeezing param- 
eter being a function of the ratio of qubit-oscillator cou- 
pling and driving strength. In turn, this state of lowest 
quasienergy has sub-Poissonian statistics. To complete 
the picture, we present the analysis of the corresponding 
dissipative dynamics for the cases when both the qubit 
and the oscillator arc damped. In presence of resonant 
multiphoton transitions, the interplay of tunneling and 
dissipation results in a stationary state which consists of 
a statistical mixture of several quasienergy states. To- 
gether with the levels involved in the respective multi- 
photon/tunneling transition, also the levels around the 
global minimum may have a large weight in the mixture. 
We show that the stationary occupation probability of 
the state with lowest quasienergy governs the lineshape 
of the nonlinear response yielding to resonant or antires- 
onant characteristics. As opposed to other anharmonic 
oscillators, the activation out of the global minimum is 
not dominantly triggered by quantum activatioii^ii, but 
instead by dissipative spin- flips (i.e., transitions between 
the two quasienergy surfaces). For Purcell-limited de- 
vices and away from a multiphoton transition, this is the 
slowest dissipative process and the corresponding rate 
can be computed by diagonalizing the Liouville operator 
numerically. Finally, we discuss the generic features of 
damped driven nonlinear oscillators in the deep quantum 
regime. 



II. COHERENT DYNAMICS IN THE STRONG 
COUPLING REGIME 

We consider a state-of-the-art circui t^^'^^'^^ or 
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cavity QED set up in the strong coupling 



regime. The cavity is driven by an external periodically 
time-dependent field. We note that this naturally im- 



plies the coupling to external modes which intrinsically 
renders the atom-cavity system an open quantum sys- 
tem. However, the strong coupling regime is character- 
ized by an atom-cavity coupling which dominates over all 
dissipative processes. In this limit, intrinsically coherent 
phenomena play a crucial role and the dissipative pro- 
cesses can be considered as small perturbation. Hence, 
in order to lay the basis for the weakly damped, dissipa- 
tive dynamics, we first focus on the coherent dynamics of 
the atom-cavity system. This paves the way to the full 
analysis of the open system which will be carried out in 
Section IIVI Hence, at this stage it is not necessary to 
restrict ourself to a specific configuration, e.g., a one- vs. 
a two-sided cavity or an circuit vs. an optical cavity. 

We model the cavity as an harmonic oscillator with 
frequency Wr which is characterized by the the ladder op- 
erators a and a' and which is coupled with strength g to 
a qubit, modeling an (artificial) atom with two only rele- 
vant quantum states, with equal resonant frequency. The 
qubit is described in terms of the Pauli operators Oj=x,y,z 
and the oscillator is modulated by a (classical) time- 
dependent field with frequency Wox and field strength /. 
The total Hamiltonian thus reads {h~l) 

H = uA a^a + -^ + -\ + g {a^ + a) Gx 

+/ (a^ + a) coswcxi • (1) 

In the frame rotating with Wcx and for the detuning 
5u) = Wr — ujcx,9,f ^ Wr, wc perform a rotating-wave 
approximation (RWA)^ and obtain the Hamiltonian of 
the driven JC model as 

H = Suj I a' a H — cr+o"_ J +— (a%_ + aa^^ — [a^ + a) , 

(2) 
with a± = CTa: ± iay. Formally, the undriven JC model 
has the quasienergies (n = 1, 2, . . . ) 



eo =: , e„,± ~ nScj ± gy/n , 



(3) 



which follow from diagonalization of the Hamiltonian in 
Eq. ([2]) for / = 0. The quasienergy states can be ex- 
pressed in the product basis of oscillator eigenstatcs {\n)} 
and qubit eigenstatcs {\g), |e)} as 



I0o) 

\(f>n+) 



|0,5), I0n-) 
1 



V2 
(|ri- l,e) + \n,g)) 



{\n- l,e) - \n,g)) 



(4) 



We will refer to the two latter as dressed n-photon or 
Fock states with two spin directions ±. The 0— photon 
quasienergy level crosses the A^-photon level at 



(5) 



Suj = ±g/VN 



(for spin =p respectively). For / 7^ 0, the crossings turn 
into avoided crossings. At such an anticrossing, the zero- 
photon state |(/)o) and the A^-photon dressed state \(J)n±) 




FIG. 1. Quasienergy spectrum (in unit of ujr) of the driven 
Jaynes-Cunimings-Hamiltonian, Eq.[2l for the parameters g — 
0.026a;r,/ = 0.004ajr. 



display Rabi oscillations with the Rabi frequency fl^ 
given by the minimal splitting of the two quasienergy 
levels. These Rabi oscillations represent the A^-photon 
transitions. In order to have well separated resonances, 
we consider here the strong coupling regime g 3> /. 
There, one would naively expect that standard pertur- 
bation theory in the driving / (and its generalization to 
quasidegenerate levels at resonance) would yield a correct 
description of all the quasienergy states which are rele- 
vant in the low energy dissipative dynamics considered 
in Section llVI However, in the remainder of this section 
we will show that for Slj of the order ^//^ standard per- 
turbation theory with respect to / fails to describe the 
lowest quasienergy part of the spectrum. In Section IIVI 
wc will show that in fact all the quasienergy states play 
an important role in the dissipative dynamics close to a 
multiphoton resonance in the stationary state. 

A first illustration of the breakdown of perturbation 
theory can be obtained by numerically diagonalizing 
the Hamiltonian in Eq. ([2]). We use the Hilbert space 
spanned by the first 100 dressed Fock states for each spin 
orientation. The result for the quasienergy spectrum for 
g = 0.026wr and / = O-OOiajr is shown in Fig. [H We 
note that wc focus on positive detuning Suj > 0, since 
the opposite case trivially follows from 5uj — > —Slu and 
|0n±) — ^ exp[— i7ra^a]|(/)„±) yielding H — > — iJ. From Eq. 
([3]) it is clear that a signature of the validity of stan- 
dard perturbation theory in the driving amplitude / is a 
constant slope of the quasienergy levels en,± as a func- 
tion of Suj. When the quasienergy splitting between two 
states is small, i.e., of the order of /, the two levels are to 
be considered quasidegenerate and an avoided crossing 
might occur. We clearly observe these patterns almost 
everywhere in the spectrum, even for levels correspond- 
ing to a large photon number that would appear in the 
region of large quasienergies not shown in Fig. [TJ The 
notable exception is the region of small quasienergy and 
small detuning. There, the slope of the levels with small 
quasienergies is not constant and no clear avoided level 
crossing appears. This clearly points to a breakdown of 



perturbation theory in / and it is easy to carry out a 
heuristic argument for this feature as foUows. 

Let us consider e„_ as a continuous function of n, 
which is zero for n = 0. It has a minimum Emin = 
~g'^/{26uj) at n,„in = g'^ / {AStu'^) and vanishes again for 
n = g'^/Suj'^. Moreover, g'^/Stu'^ can be regarded as the 
number of states in the quasicnergy interval emin < e < 0. 
The average quasicnergy spacing e in this interval is thus 
Suj/2. Moreover it is clear that the levels tend to accu- 
mulate close to the minimum and that e overestimates 
the local level spacing there. The matrix elements of the 
driving term in this region are of the order ^Jn-axinf . We 
can conclude that the driving term is no longer a per- 
turbation in the lowest energy part of the spectrum if 
V^-min/ ~ e or (5a; ~ y/7g- This rough estimate is con- 
firmed by the numerical data. In fact, for / = 0.004a;r 
and g = 0.026a;r as in Fig. [1] we find ^/Jg « 0.01, which 
correspond approximately to the region where the slope 
of the levels in not constant. 



III. QUASIENERGY LANDSCAPE AND 
DYNAMICAL BISTABILITY 



Motivated by the failure of ordinary perturbation the- 
ory, we next formulate a different pcrturbative approach. 
For this, it is convenient to switch to the picture of 
dressed qubit states, formally achieved by the unitary 
transformational 



R = exp 



-Svr 



S^a^a + (7+cr_/4 



(aV- 



acr+J 



(6) 



It maps the JC eigenstatcs Eq. (JH) into product states 
according to 



|0,.9)^|0,g). 



W.g), \(t)n+) 



l,e). 



The purpose of this transformation is that the undriven 
JC Hamiltonian becomes diagonal in the qubit Hibert 
space and assumes the form 



H = \5uj\ ( a* a + -a+a^ 



gf^z 



ata+ ^CT+cr- , (7) 



while the ladder operators remain unaffected to lowest 
order in the photon number, i.e.. 



WaR = a + 0{n-'^^) 



1/2N 



(8) 



The expression for a follows by expanding the matrix 
elements 

{ng\a\mg) = {ng\R'' aR\mg) = ((^„_|a|(^„i_) 

= /+(")'5)i+lm , 

{ng\a\me) = {ng\R^ aR\me) ~ {(j)n-\a\(j)m+) 

= f-{n)5„+i,n , (9) 



with /±(n) = (^/?^ ± Vrr+T)/2. Note that /_(n) is 
typically small, e.g., for n = 1 (it is the largest case). 



az = 1 



az = -1 



1.2 




-0.35 



FIG. 2. 3D plot of the quasienergy surfaces Q{X, V) for P > 
and f/g = 0.154 (Note that this yields an effective tempera- 
ture of Tofi = 0.65, see text). The quasienergy surface for the 
dressed spin orientation e (ctz — +1) is a monostable reversed 
cone. The one for the opposite dressed spin g (a^ = —1) is 
bistable and is shown as a density plot as well. In the den- 
sity plot, the oriented black solid line indicates the separatrix 
between the two domains of attraction. 



/_(1) sa —0.2 and it approaches zero as — 1/(4y^). On 
the other hand, f+{n) ^ -y/n + O^l/^/n). This also il- 
lustrates that only higher order terms depend on spin 
flipping operators. Hence, it follows that both the driv- 
ing and the coupling to the bath induce spin flips, but 
only as higher order processes. 

Next, we introduce the rescaled rotating quadrature 



X 



Sy 



+ a). 



r^i 



/|(«^ - a) 



(10) 



with A = \duj\'^ / g"^ . Note that the terms neglected in 
Eq. ([8|) arc of higher order in A. By plugging the ro- 
tating quadrature into the transformed Hamiltonian Eq. 
([7]) and neglecting all higher order terms, we obtain 
H c^5ujX-^Q{X,V) with 



QiX,V) 



X^ 
2 



2 



X^ 



'El 

2 



/ 
V2.g 



X. (11) 



Eventually, rescaling the time by r = 
effective Schrodinger equation 



l/(5w, yields the 



(12) 



with the commutator [X ^V] = i\ (which follows from 
the definition in Eq. (jlOp ). Thus, we can interpret A 
as a rescaled effective Planck constant, X and V as 
canonically conjugated operators and Q{X,V) as two 
quasienergy surfaces in phase space for the two opposite 
dressed spin orientations. They are visualized in Fig. [21 
The quasicnergy surface for the dressed spin g {az — —1) 



has the overall shape of a Mexican hat. The drive in- 
duees a finite tilt of the surface, generating a saddle point 
at T'sad = and Xgad = (1 ~ f/g)/V^ and a separa- 
trix, shown in the density plot in Fig. [2] as the oriented 
black solid line. It divides the surface into three domains: 
i) a potential well around the quasienergy minimum at 

^min = and A'niin = ~(1 + //.9)/v2, fi) an internal 
dome around the inner maximum at Af = 0, 7^ = and, 
iii) an external surface. For quasienergies lying above 
the saddle point and below the maximum, there are or- 
bits coexisting on the domains ii) and iii) and thereby 
define a dynamical bistability. The surface for the oppo- 
site dressed spin orientation a^ = -1-1 is a less interest- 
ing monotonous function and is shown in the 3D plot in 
Fig.H 

When the motion is quantized, the (quasi-) energy lev- 
els become discrete. For small A, we expect the dynam- 
ical behaviour to be semiclassical and we can associate 
each quasienergy level to an allowed orbit. In principle, 
a full WKB treatment is possible but it goes beyond the 
scope of this paper. In order to illustrate the semiclas- 
sical features of the driven Jaynes-Cummings model, we 
show the quasienergy levels (rescaled by X/5uj) as a func- 
tion of A in Fig. |31 They are obtained by numerically 
diagonalizing the complete Jaynes-Cummings Hamilto- 
nian. 

We can associate the levels below the quasienergy 
saddle point Q{Xsad,0) = —(1 — f/g)'^/A (indicated in 
Fig. [3] as a green dotted line) to orbits localized in the 
quasipotential well. Notice that close to the minimum 
Q{Xmin, 0) = — (l+//(7)^/4, the orbits are approximately 
equally spaced and they become denser while approach- 
ing the saddle point. In fact, semiclassically the level 
spacing is given by Xuj{q), where uj{q) is the classical or- 
bit frequency for the (rescaled) quasienergy q. Close to 
the quasienergy minimum the potential is harmonic and 
Lu{q) varies slowly while close to the saddle point, it tends 
quickly to zero. 

For energies above the saddle point and below the max- 
imum, we can associate the quasienergy levels with orbits 
in the internal dome and in the external surface. The or- 
bit with vanishing quasienergy corresponds to an orbit 
close to the quasienergy maximum. For certain discrete 
values of A, an orbit on the external surface might have 
vanishing quasienergy as well. Hence, we can interpret 
the multiphoton transitions, which occur at the resulting 
avoided level crossings, as tunneling transitions between 
two states in the internal dome and the external well, re- 
spectively. In the limit of g ^ /, we can read off the res- 
onant tunneling condition from the perturbative result in 
Eq. ([S|), yielding the condition for iV-photon transitions 
to be A = 1/7V. More tunneling transitions for negative 
quasienergy are indicated by the orange circle in Fig. [31 

Above the minimum no avoided level crossing is 
present, several exact crossings (not shown) appear in- 
stead. The corresponding orbits have opposite dressed 
spin and do not show level repulsion because the two 
quasienergy surfaces are effectively decoupled. 



to 




FIG. 3. Rescaled quasienergy spectrum of the driven Jaynes- 
Cummings Hamiltonian in Eq. ^ for the parameters f/g — 
0.154 (solid black lines). The green dotted line indicates the 
saddle point quasienergy. The red dashed lines correspond 
to the approximate {X/6oj)En, see text. The blue dotted- 
dashed line marks the rescaled quasienergy £o of |0(/)). The 
orange circle highlights avoided quasienergy level crossings 
(multiphoton transitions) for large photon numbers A''. 



A. Quasienergy states at the well bottom 

Next, we study the spectrum close to the bottom of 
the quasienergy well quantitatively. In fact, as it turns 
out, the states localized in this region play an impor- 
tant role in the stationary dissipative dynamics close to 
a multiphoton transition. 

The most simple approach consists in expanding the 
quasienergy in Eq. (jlip close to its minimum. In this re- 
gion, the quasienergy surface is to lowest order harmonic 
and follows 



Q{X,V) = Q{X^in,V, 



1 



mini ' mmf 



2mcff 



iv-r„ 



--moffW* {X - Xnunf 



(13) 



with effective mass rrics = (/ + g)/f and with fre- 
quency u!* = \/ f /{f + g)- When the zero point energy 
Xlo* /2 is much smaller than the quasienergy well depth 

AQ EE Q{Xsad,'Psad) ^ Q{Xmin,'Pmin) = f/g, Only a feW 

quasienergy states are localized in the well. The har- 
monic expansion yields the quasienergies 



E„ 



5ujX 



QiXminiO) + AtJ* ( n + 



(14) 



which are determined up to ©(A^)^. They are shown in 
Fig. [3] as red dashed lines and they almost coincide with 
the exact results for small A. 

Up to leading order in A, the corresponding wave func- 
tions are given by 



m 



R-^—=D{X„ 
Vnl 



.)Sir*)\Og). 



(15) 



They are obtained by applying to the vacuum: i) the 
squeezing operator S{r*) = exp[r*{a'^ — a^'^)/2\ with 
squeeze factor 



ln[mcffW*] ln[l+.g//] 



(16) 



ii) the translation D{X^\n) ~ cxp [iVXynm/M to the mini- 
mum, iii) the creation operator 6^ — a^ cosh r + a sinh r — 
e^Xynin/V2X, and iv) R~^ to return to the bare qubit 
picture. With this, one can compute the Fock-state 
representation^^ of the quasicncrgy states and all expec- 
tation values at leading order. 



B. Quantum squeezed state and sub-Poissonian 

statistics 

Of particular interest is the state jV'o) at the bottom 
of the quasienergy well. It is defined, when the zero 
point energy Xui* /2 is smaller than the quasienergy depth 
AQ = f/g, corresponding to the region Slo < g\/f/g 
(= 0.016 for the parameters used in Fig. [1]). Note that 
the analytical result for Eq in Eq. (|14p almost coincides 
with the exact one in this region. Most importantly, \^Pq) 
exists for any finite driving (providing that the detuning 
is small enough), and is clearly absent in the undriven 
case. In fact, for / = 0, the tilt of the quasienergy sur- 
face vanishes and no quasienergy well is developed. This 
illustrates further that simple perturbation theory fails 
for arbitrary small driving close to zero detuning and im- 
proves the estimate of the critical detuning below which 
it breaks down. 

For weak but finite driving, the well is still very shallow 
in the momentum direction, allowing for large momen- 
tum fluctuations of a state conflned in it. For increasing 
driving, the well becomes deeper and more symmetric. 
Hence,|V'o) i^ amplitude squeezed and the squeezing de- 
creases for increasing driving, as it can also be read oflF 
from Eqs. ([75]) and (fTHl). 

As a consequence, jV'o) has also sub-Poissonian statis- 
tics and shows photon antibunching. This can be seen 
by computing its average photon number fi* = (■(/'ol"-|'/'o) 
and its variance (An*)^ = {tpQ\{n — n*)'^\iJjQ) by means of 
Eq. (dH). We find 



X^,^_{.f + 9f 



and 



2A 45^2 ' 



*\2 -2r*-* -2r* (j + 5) 



(An ) ~ e n = e 



ASlu^ 



(17) 



(18) 



Thus, in the semiclassical limit A — > 0, the 2-photon cor- 
relation function becomes 

.^^nO)^l + ^^^^^^<l, (19) 



Note that the above formulas can be regarded only as 
leading-order estimates, because they are of order 1/A 
and other finite contributions are not taken into account. 

Here, we did not include the Fock representation of 
the semiclassical states because it is rather cumbersome 
and scarcely illuminating for the purpose of this work. 
However, it is important to keep in mind that the semi- 
classical states are a superposition of dressed Fock states 
with a given dressed spin orientation and photon number 
n of the order of n, given in Eq. (|17p. For the discussion 
of the dissipative dynamics below, we remark that close 
to a A^-photon transition, all the dressed Fock states that 
contribute significantly to the superposition correspond 
a photon number n < N with negative quasicnergies, 

£n<N- < 0. 

Finally, we note that the Golden Rule estimate for the 
rate for coherent spin flip transitions yields zero, because 
no overlap between the quasienergy states localized in 
the well (which has negative quasienergy £"„ < 0) and 
that with opposite spin (which has positive quasienergy 
£„-(- > 0) exists. 

In this work we consider identical resonant frequen- 
cies for the qubit and the cavity. However close to 
the quasienergy minimum the effect of a finite detun- 



ing (5 



: UJn 



- ojc of the resonant frequencies Wq and Uc for 



the qubit and the cavity might be negligible even if the 
circuit cavity QED is in the dispersive regime when few 
photons are present. The detuning S induce an additional 
term S'^X/{4:g'^) inside the square root in Eq. (fTT|) . This 
further contibution is negligible close to the well bottom, 
if 5^X/g^ < 2A'^i„ = (1 + f/gf- Note that this condi- 
tion can be fuUfiUcd even li 5 ^ f,g a.s in Rcfs.— 1^. In 
this references, the state \^q) has been denoted as bright 
state. 



C. Dark state and multiphoton transitions 

Another state playing a key role in the dissipative dy- 
namics is the state obtained by starting in the JC ground- 
state |0, g) and by adiabatically switching on the driving. 
This is the state with zero quasienergy and is shown as 
blue dotted-dashed line in Fig.|31 We denote it as |0(/)). 
It is naturally favored by photon leaking, the most impor- 
tant dissipation channel in many set-ups. For weak driv- 
ing, it has vanishing average photon number and there- 
fore it is not accurately described by our semiclassical 
approach. 

This state has been investigated in Ref. [30 for 5uj ~ 0. 
It is a squeezed state which tends to follow the driving 
by rotating its spin by the angle 6 = arcsin f/g around 
the the y— axis and squeezing its amplitude fluctuations 
with squeezing factor tq = -[ln(l - p/g^)]/4^. Note 
that, as opposed to r* (the squeezing factor of |'0*)), 
To increases for increasing driving. In addition, it has 
supcr-Poissonian statistics and shows photon bunching. 



In fact, 



no = (0(/)|n|0(/))/=sinhVo, 



(20) 



{Anof={0{mn^nof\0{f)}^2cosh^rosmh^ri)21) 
and the 2-photon correlation function is 



ff('HO)^l + 



(An) 



> 1, 



(22) 



Since it is not possible to diagonalize analytically the 
driven JC Hamiltonian for Suj — 0, a systematic pertur- 
bation theory in Suj is not possible. Fortunately, for weak 
driving, we can rely on perturbation theory in the driving 
/ to study this state. Away from a multiphoton transi- 
tion, we can compute the transmission (0(/)|a|0(/)) by 
ordinary perturbation theory in /, yielding 



(0(/)|«|0(/)) 



1 



1 



/ 

A ySio + g ' 5lo — g 



(23) 



Hence, the oscillation of the transmission is in phase with 
the driving for 6u) < g. 

At a multiphoton transition, i.e, for A ~ 1/N , the cor- 
responding multiphoton state is a superposition of |0(/)) 
and a semiclassical state \N{f)), which exists on the ex- 
ternal part of the quasienergy surface. This state is ob- 
tained by starting from a dressed state \(t)N+) (l^w-) 
for 6lu > 0) and switching on adiabatically the driving. 
The states |0(/)) and \N{f)) display Rabi oscillations. 
The corresponding Rabi frequencies can be computed by 
means of Van Vleck perturbation theory as 



fii 



^/2' 



n. 



%/2. 



P 



n.^ 



3^/3 

24 g2 



(24) 



We do not give the general formula because it is rather 
cumbersome and scarcely illuminating. In fact, it is de- 
rived assuming that perturbation theory is valid every- 
where in the spectrum, which is generally not true, not 
even for f -^ g. Our approach does not allow a fully- 
fledged semiclassical calculation of the Rabi frequency. 
Nevertheless, one main advantage is that it describes very 
elegantly how the physical quantities are rescaled while 
the effective Planck constant changes. With this, we pos- 
tulate that the Rabi frequency follows as 



n^ ^ exp [-S{f/g)/X] = exp [-NS{f/g)] 



(25) 



with S{f / g) being an unknown function of f / g and 
indicating logarithmic precision. 



the timescale Tr — I-k jbJr-, governing the system dynam- 
ics. For instance, for the transmon architecture reported 
in Ref. [25I T\,2 ^ Ia's and T^ ^ 1 ns. For optical cavities, 
quality factors of Q ~ 10^° are possibleS^. When deco- 
herence and dissipation are induced by electromagnetic 
environmental fluctuations with a smooth spectral den- 
sity and when all the time scales governing the different 
dissipative processes exceed typical bath-intrinsic corre- 
lation times, a Markovian dynamics is expected. The 
simplest Markovian master equation (MME) for the re- 
duced density operator pit) of the qubit-plus-oscillator 
system is of Lindblad form and incorporates oscillator re- 
laxation, qubit dephasing and qubit relaxation. Its stan- 
dard form is given by p = £ [p] , with the Liouvillian 



L{p\ = -i\H,p\^KD\a\ 



7i2?[a-/2] + ^I?[<T,],(26) 



where 'D\0\ is the Lindblad damping superoperator 
V{0\ = ([Op,0^]-H[0,pOt])/2. The three damping terms 
describe: i) the photon leaking out of the oscillator at 
rate K, ii) intrinsic qubit relaxation at rate 71, and iii) 
pure qubit dephasing at rate 7,^. This phenomenological 
master equation follows by modelling the environment as 
three independent harmonic baths, each held at thermal 
equilibrium at the same temperature T, provided that 
hK,h'yi,h'y^ <C k^T,fiWr- In Eq- (|26p . we have implic- 
itly assumed a one-sided cavity with a single input and 
a single output port, yielding a single dissipative channel 
for the cavity. For a detailed description of a nanocircuit 
and an optical set-up implementing this model, we refer 
the reader to Ref. [2^ and Ref. ,23, respectively. 

Moreover, one has to assume fcsT <C ftw^- This con- 
dition together with the ordinary assumptions for the 
RWA {Sijj,g,f <C uJr) implies that the environment acts 
as a perfect energy sink. However, quasienergy and not 
energy is the good quantum number for driven systems. 
Since the quasienergy is defined in a rotating frame, en- 
ergy emission in a static frame (energy leaking into the 
environment) can appear as quasienergy absorption in 
the rotating frame. This leads to counterintuitive effects, 
such as the Unruh effect for a constantly accelerated rel- 
ativistic system^ and quantum activation for the quan- 
tum Duffing oscillator—. In the following, we will illus- 
trate how our approach in terms of the quasienergy sur- 
face provides an intuitive physical insight into the dissipa- 
tive dynamics even when many different quantum states 
are involved. 



Potential well escape: quantum activation vs. 
spin flips 



IV. DISSIPATIVE DYNAMICS 



State-of-the-art nanocircuit QED set-up» 



15.16.18.25 



are 



characterized by weak damping, implying large quality 
factors of the order of Q ~ 10^ — 10^ and qubit dephasing 
and relaxation times (Tj* and T2) being large compared to 



Away from any multiphoton transition and for g ^ f, 
|0(/)) is characterized by a vanishing average photon 
number. If the system acts as an energy sink, this state 
has an infinitely large lifetime because a photon cannot 
be emitted neither from the system nor from the environ- 
ment. Hence, |0(/)) will be dominantly populated in the 



stationary state. Conversely, at a multiphoton resonance, 
the system can escape from the 0— photon state |0(/)) 
via resonant tunnehng to the dressed A^— photon state 
|iV(/)). Subsequent relaxation causes the energy to leak 
out from the system via transitions to states with lower 
quasienergies. This leads to the occupation of all the 
quasienergy states with negative quasienergies. In fact, 
these states are a superposition of dressed Fock states 
with photon number n < iV, as detailed in the previous 
section. Hence, when Scj < g y/f/g, energy leaking leads 
to the occupation of those quasienergy states which are 
localized in the quasienergy well. 

From the basin of attraction of the quasienergy mini- 
mum, the system can escape either i) by climbing up the 
quasienergy well by quantum activation, or ii) by means 
of a spin flip transition. In fact, only one of the two 
dressed spin states is confined. Hence, after a spin flip 
transition the system can quickly decay to the dark state. 
The rate Ta for quantum activation is suppressed expo- 
nentially with the number Ny^ of states in the quasienergy 
well, following F^ ex Ke~^^'" . Note that the constant c 
does depend on the precise shape of the quasienergy sur- 
face, but not on the rescaled Planck constant A, whereas 
N^ obviously depends on A (which governs the level spac- 
ing) according to N^ ex 1/A. Hence, spin flip transitions 
are the dominant escape channel when A is small. 

The stationary state is a statistical mixture of all those 
quasienergy states with negative quasienergy, if the over- 
all rate of escape from the quasipotential well is of the 
same order as the rate of escape from the state |0(/)). 
The latter is discussed in the next subsection. 



B. Coherent vs. incoherent resonant dynamical 
tunneling 

Driving may also induce resonant transitions between 
the qubit dressed states. To be more precise, we have 
to distinguish between the regimes of (A) coherent and 
(B) incoherent resonant dynamical tunneling. The bor- 
derline between these two regimes is determined by the 
ratio rijv/rAT, where VLff is the Rabi frequency of the 
TV— photon transition and F^v is the inverse lifetime of 
the dressed A^— photon state. The latter is of the order 
of the maximum of TVk, 71, 7^. 

(A) Coherent resonant dynamical tunneling is char- 
acterized by the decay rate of the dressed iV— photon 
state being smaller than the Rabi frequency, F^r < VLm- 
Then, the system coherently tunnels back and forth with 
the Rabi frequency 0,^ several times before it signifi- 
cantly relaxes with decay rate FAr/2. All the subsequent 
dissipative transitions occur on a time scale similar to 
Tm- Therefore, the stationary solution will be a statisti- 
cal mixture of all the states with negative quasienergies. 
This is qualitatively different from the situation away 
from resonance where only the 0-photon state is signifi- 
cantly populated. 



(B) Incoherent resonant dynamical tunneling occurs 
when the decay rate of the A^— photon dressed state 
exceeds the Rabi frequency, Fjv > ^n- Then, the 
quasienergy level broadening is larger than the coher- 
ent splitting of the two quasienergy levels. The dressed 
A^— photon state strongly fluctuates on a quasienergy 
range F^v and with it the quasienergy difference to the 
0— photon state. Tunneling occurs only for those split- 
tings which do not exceed the Rabi frequency ^n- This 
occurs with probability Jljv/rAf- When the two levels are 
quasi-degenerate, a tunneling event then happens with 
probability fi^r. For incoherent resonant dynamical tun- 
neling, the total tunneling probabilty is just the prod- 
uct of the two. Hence, the system can escape from the 
0— photon state jO^) with a small total rate JI^/Ftv- In 
the incoherent regime, again two situations can arise: 

(Bl) When no state is localized in the quasienergy well 
(which roughly occurs for Su > g\/f/g, see above), i.e., 
when perturbation theory in the driving / is valid, the 
lifetimes of all the states visited during the relaxation 
transition are of the same order, namely of the order of 



F^ . This is given by the maximum of k' 



-^7^^ 



y. 



This, however, is in any case much smaller than the life- 
time of the 0-photon state, given by Fat/O^^, (remember 
that Fjv is suppressed exponentially with A'^). Thus, the 
0-photon state is maximally populated. 

(B2) In the second case, when several states are 
trapped in the quasienergy well (which occurs for Suj <C 
9yf/g)-! i-6-; when the semiclassical description is neces- 
sary, the rate for quantum activation F^ is exponentially 
small and can be of the order of ri^/Fjv. If also the 
spin flip rate is small, the state at the well bottom can 
become metastable. In this case, it is therefore possible 
that an overall small rate of incoherent resonant dynam- 
ical tunneling leads to a dramatic change in the station- 
ary distribution as compared to the situation away from 
resonance, resembling the situation of coherent resonant 
dynamical tunneling. Note that, this will be the case only 
for the first few resonances. In fact the resonant escape 
from the dark state is suppressed exponentially with A'^, 
which is not the case for the spin fiip escape from the 
bright state. 



V. NONLINEAR RESPONSE IN THE STRONG 
COUPLING REGIME 

Having discussed all required ingredients, we can next 
turn to the nonlinear response of the driven JC model, 
which can be characterized by the steady-state expecta- 
tion value of the intracavity field operator related to 



tT{ga) = Ae}' 



7 ^ Qal3a.j3a 



(27) 



The indices a, /3 refer to the basis of eigenstates of the 
driven JC Hamiltonian in Eq. ([5]). In a one-sided cavity, 
its modulus A is related to the transmitted amplitude 




co„/co. 



FIG. 4. Nonlinear response of the driven JC model: am- 
plitude A as a function of the driving frequency Ucx for 
/ = 0.0353 (a), / = O.O7I3 (b) and / = 0.17g (c). More- 
over, K = 6.1 X 10~^q, 71 = 1.4 X 10"^a, 7<4 = 1-4 X 10"*' 



g = 0.007ajr. A realistic value for ajr/(27r) is 7 GH: 



,25.38 



9, 



Atr oc A and intensity /t,- oc A^ of the input signal in a 
heterodyne measurement set-up2^>^. 

These quantities are directly accessible in the exper- 
iment. In the rotating frame, (a) < corresponds to 
an oscillation out of phase ((p = tt) with respect to the 
drive, while the in-phase oscillation ((/) = 0) is associated 
to (a) > 0. Note that 6uj — >■ —5u) and Qap — >■ g*fl, since 
the unitary transformation exp [— iTra^a] yields C -^ C* . 
Moreover, apa -^ —o-pa and all matrix elements are real. 
We can conclude that the nonlinear response for Suj < 
follows trivially from the one for Slj > according to 
(a) — ?> —(a)*, i.e. A ^ A and ip ^ n — (p. In addition, 
due to the RWA, the master equation can be written in 
terms of the ratios f/g, n/ g, 7i/.g , 70/3 and A only. 
Thus, also (a) depends only on these quantities. 

A straightforward numerical solution of the stationary 
limit pit) = allows to numerically calculate the modu- 
lus A as a function of the external modulation frequency 
Wcx- The result for experimentally realistic parametera^^ 
is shown in Fig. |4l For weak modulation, two large fun- 
damental (anti)resonances appear symmetrically with re- 
spect to Wcx = Wr, which mark the supersplitting of the 
vacuum Rabi resonance (in Fig. 2^', we show only the 
regime of Wex < Wj). These largest resonances (which are 
in fact antiresonances) are associated to the 1-photon 
transitions. They can be describe d^ ^'^^ by an effective 
two-level model involving the 0- and the 1-photon state 
only. At resonance, both states are equally populated 
and oscillate with opposite phase. This overall yields zero 
response exactly at resonance. Slightly away from reso- 
nance, one of the two states is slightly more populated 
and a finite response arises. Further away from the res- 
onance, the response again approaches zero as discussed 
above (see below for a quantitative evaluation). Hence, 
overall, the lineshape of an antiresonance arises. The 1- 
photon transitions are special situations since only two 
quasicnergy states are involved and no contributions of 
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FIG. 5. Solid lines: Nonlinear response of resonator in the 
circuit QED setup based on the transmon architecture as re- 
alized in Ref.m, i.e., uj/{2-k) = 6.920 GHz, g/n = 94.4 MHz, 
K/(27r) = 300 kHz, //(27r) = 8.304 MHz, 7i/(27r) = 98 kHz 
and 70 = 0. Dashed lines: same parameters, but without di- 
rect qubit relaxation (71 = 0). Shown are (a) the amplitude 
A, (b) the phase cj>, (c) the quasienergy spectrum, and (d) the 
population P* of the lowest quasienergy state IV'o)- The red 
dot-dashed line in (a) shows the result for the population of 
the 0-photon state being assumed as close to one, see text. 



other quasisenergy states with smaller photon numbers 
exist. Increasing the modulation strength, further reso- 
nances and antiresonances appear, see Fig. 2]b) and c). 
The antiresonances turn into resonances when the driv- 
ing is further increased. 

In order to illustrate the underlying mechanism, we 
show another example of the modulus A (Fig. [S^) and the 
phase <f (Fig.[5j)) of the oscillator response, together with 
the associated quasienergy spectrum (Fig. [5t) and the 
occupation probability P* of the state \'(p^) with lowest 
quasienergy (Fig. [St). For the solid lines in Fig. [SJ we 
have chosen the same parameters as in the fits of Ref. 
I25I . As it turns out, in this experiment, the qubit pure 
dephasing is negligible (7^ = 0). We have also plotted 
(dashed lines) the case of pure Purcell dissipation (no 
intrinsic qubit relaxation 70 = 0, 71 = 0, only finite 
resonator damping). 

First of all, we confirm that the resonances and antires- 
onances occur in correspondence with the avoided cross- 
ings in the spectrum. As already discussed in the previ- 
ous section, in absence of resonant tunneling, i.e., away 
from any resonance, the occupation probability of the 
0-photon state is close to one (under the conditions dis- 
cussed above), i.e., poo ~ 1. Then, from Eq. (^5)) follows 
that (a) w aoo = (0(/)|a|0(/)), which is shown as red 
dashed line in Fig.[St). The first antiresonance occurs for 
Wox ^ ^r ^ 9 and has already been discussed above. Figs. 
[5^) also shows the 2— ,3— ,4— photon resonances and the 
5— ,6— photon antiresonances. The overall phase of the 
stationary oscillations changes in presence of a resonance 
(Figs. [5)3), which points to a significant population of a 
state oscillating out of phase. In fact, the occupation 
probability P* of the state 1-00) with lowest quasienergy 
displays peaks at such resonant frequencies, see Fig. [5Jl). 
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This phenomenology is consistent with the picture of 
the dissipative dynamics drawn in the previous section. 
At a multiphoton resonance, in presence of coherent reso- 
nant dynamical tunneling or when the rate of incoherent 
resonant dynamical tunneling is of the order of the rate 
for an escape from the quasienergy well, a sizeable occu- 
pation probability of the states confined in the quasipo- 
tcntial well results. Those quasienergy states oscillate out 
of phase and with a large amplitude, yielding peaks in the 
nonlinear response and an overall out-of-phase response 
characteristics. In the opposite limit of incoherent reso- 
nant dynamical tunneling and when no deep quasienergy 
well exists (i.e., in the perturbative regime) or when the 
escape rate is smaller than the incoherent tunneling rate, 
the 0-photon state is again dominantly populated and 
a small but finite occupation probability of at least one 
state oscillating out of phase leads to a reduction of the 
resonant nonlinear response and thus to a dip in the line- 
shape, but does not necessarily change the phase. We 
note that the underlying mechanism is exactly the same 
as for the multiphoton transitions in the quantum Duff- 
ing oscillator—"—. 



VI. PURCELL-LIMITED SET UPS 

The comparison of the solid and dashed lines in Figs. 
[S^, b, d) shows that a weak intrinsic coupling of the 
qubit to the environment docs not change the underlying 
physics qualitatively. In particular, transmon qubits are 
Purcell-limited in the resonant regime^. It is therefore 
interesting to consider the special case of pure Purcell 
dissipation 71 = 70 = 0. In this limit, the dissipative dy- 
namics is governed by the simplified Liouville operator 

C[p] EE -z[H,p] + |([ap,at] + [a, pat]). (28) 

We distinguish between two different regimes: i) when 
only few photons are exchanged and no state is localized 
in the quasienergy well, a perturbative analysis (in the 
driving) is in order, ii) In the opposite case, our physical 
picture in terms of quasienergy surfaces is necessary to 
discuss the dissipative dynamics. 



A. Small photon number: perturbative regime 

When no state is localized in the quasienergy well and 
away from a any multiphoton resonance, the quasienergy 
levels match the unperturbed result of Eq. ([3]) and the 
corresponding states can be identified with the dressed 
states \(j)n-)- 

Close to the A^— photon resonance, two scenarios are 
possible, as discussed above: (A) coherent resonant dy- 
namical tunneling (i.e., when the lifetime Tn of the 
A^— photon state \4>n±) is much smaller than the Rabi 
frequency Qat), and, (B) incoherent resonant dynamical 
tunneling (i.e., when T^ ^ ^n)- 



When the dynamical tunneling is coherent (F^r <C 
^n), the system tunnels many times between the states 
\(j)o) and \(j)N-) with period 2-k/^n before it substantially 
decays to \4>m-i-)- The rate of this decay can be ob- 
tained in secular approximation as FAr/2 = k{N—1/2)/2. 
Subsequent decays from \4>n-) to |0„_i_) occur along 
the ladder A-l^A^-2^... ^1^0 with de- 
cay rates 2?n-,(,i-i)- = (v^ + y/n — 1)^k/4 for n 7^ 1 
and I'l-^o = k/2. Hence, the rate from the 1-photon to 
the 0-photon state is smallest. Note that the probabil- 
ity of a decay to a state with opposite dressed spin is 



small, i.e., 2?, 



(n-l)-,n+ 



/2?(„-l 



l/[16n(n-l/2)] 



Thus, in this scenario, the occupation probability of pn 
is the largest. It can be computed (once the spin-flips 
are neglected) by straightforwardly solving the master 
equation, which yields 



N-l r-2 



Pi 



,= [l + f-^iN)+J2i 



in) 



(29) 



where /+ is defined below Eq. (j9]). Therefore, the contri- 
bution /Oil ail with 



an 



/ 3 + 2^/2 



/ 



2(ei_-£2- 



4(ei_ -£o) 



(30) 



is the largest in Eq. ((27|) . For all multiphoton transitions 
with photon number A < 6, ei_ < e2-,£o, resulting in 
the contribution pnan < and thus in an overall sta- 
tionary oscillation which is out of phase with the modu- 
lation. 

In the example considered in Fig. [!](/ = 0.057(7, 
K — 0.0009(7 and g = O.OOToj,.), this scenario is realized 
close to the 2— photon transition. The ratio between the 
2-photon Rabi frequency ^2 = \/2p/g (see Eq. ((24)) ) 
and the lifetime of the 2— photon dressed state r2 = 3/2k 
is l^s/Fz = {2V2/3) f/(gK) = 3.4 and the tunneling is 
coherent. As expected, the overall behavior changes from 
an in-phase oscillation with amplitude |aoo| (red dashed 
line in Fig. IH]) to a larger out-of-phase oscillation of the 
order of janj (blue dotted line). Moreover, the occu- 
pation probability pn displays a peak (Fig. HJl), whose 
height is close to the approximate value computed in Eq. 
([29]). It is represented as a blue star in Fig. [6jl). 

In the opposite limit of incoherent resonant dynamical 
tunneling {Tpj 3> ^In), the population of \(j)o) is ~ 1, be- 
cause relaxation to the 0— photon state is more efficient 
than any escape from there via tunneling. Dissipation 
then completely washes out the resonance, and the re- 
sponse is identical to that away from resonance and thus 
is in phase with the drive. 

In the intermediate regime when F^r ~ ^n, a small 
population pn emerges, contributing pnan with an < 
given in Eq. (PU)). which is negative for A^ < 6 (in 
the perturbative regime). This also leads to a reduced 
response with an antiresonance. This scenario is realized 
close to the 3— photon transition in Fig. [S]). In this case, 
the Rabi frequency is fl^ = {3'^/^/2'^)f/g^ and T3 = 
(5/2)k, so that fls/Ts = 0.24. 
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FIG. 6. Nonlinear response of the driven JC model: (a) am- 
plitude, (b) phase, (c) quasienergies, and (d) population of 
the lowest quasienergy state for / = 0.057(?, k = 0.0009(? and 
g — 0.026ajr. Dashed red (dotted blue) line in (a,b,c): lowest- 
order result for |0(/)) (|l(,f)>). 



FIG. 7. Nonlinear response of the driven JC model: (a) am- 
plitude, (b) phase, (c) quasienergies, and (d) population of 
the lowest quasienergy states \tp^) [n = 0,1,2) as a func- 
tion of the rescaled Planck constant A for the parameters 
f/g = 0.171, «:/.g = 6.110"-'ajr, 71,7^ = (PurceU limit). 
Green dotted line in (c): quasienergy of the saddle point. 



B. Large photon number: semiclassical regime 

When several photons are contained in the resonator, 
a picture in terms of the dissipative semiclassical dy- 
namics emerges. As will be shown below, a separation 
of relaxation time scales exists, which separates fast in- 
trawell from slow interwell relaxation. We first focus 
on the intrawell dynamics. Close to the minimum, the 
wave functions are given by Eq. (jlSp . We can compute 
the corresponding dissipative transition rates by plugging 
them into Eq. (j28p . In this limit, dissipative transitions 
occur only between nearest neighbors, with the rates 
'Dn-i.n = uncosh r and I?„^„_i ~ Knsinh r. Here, the 
detailed balance condition is fulfilled. Hence, when the 
system is initially in a state with a large photon number 
n, it falls with large probability into the basin of attrac- 
tion of the quasienergy minimum. This process consti- 
tutes intrawell relaxation. When furthermore damping is 
smaller than the intrawell level spacing, i.e., n <?C g\uj* 
(see Eq. (HH)), detailed balance is retained and deter- 
mines an effective Boltzmann distribution 



p: = P* 



-"^o: 



(31) 



with the effective inverse temperature jS^a = 2 In coth r* 
being defined in terms of the squeezing parameter r* . We 
emphasize that this link between effective temperature 
and squeezing can be generalized to any driven quantum 
system with a smooth quasienergy surface and coupled 
linearly to a bath (e.g., the linearly and the paramet- 
rically driven Duffing oscillator). It can be easily gen- 
eralized to finite real temperatures T > as well. It 
turns out that the zero temperature limit applies when 
sinh r is much larger than the bosonic occupation num- 
ber h(ujcx/T) of the bath taken at frequency Wex- In the 
opposite limit, /3eff = ^ex/T. Since we include here only 
photon leaking, i.e., Wox 3> T, the effective temperature 
is still small. 



On the large time scale, the system can escape from the 
basin of attraction of the quasienergy minimum with a 
small forward rate k^. Overall there are two mechanism 
contributing to the escape: i) The system can climb up 
the quasienergy well by quantum activation. As already 
discussed in Section lTVl the corresponding rate Ta is sup- 
pressed exponentially, following Ta c>c Ke~'^'^. Moreover, 
ii) the bath induces spin-flip transitions with rate Fg even 
in absence of an intrinsic spin-bath coupling. This is a 
consequence of the dressing of the spin. From Eqs. (|8]) 
and (pS)) . it follows that Tg ex. A. Therefore, this mech- 
anism is dominant for small A and the rate fc+ vanishes 
according to A:+ oc A. In real set-ups, a lower bound to 
k^ is given by the small intrinsic couplings 71 and 7^ of 
the spin to the bath. 

Once the system has left the basin of attraction of 
the bright state \iPq), it quickly decays to the dark state 
|0(/)). From there, it can escape with a backward rate 
k- on a large time scale. The stationary populations 
of the intrawell states, which oscillate out of phase, and 
the 0— photon state, which oscillates in phase with the 
modulation, are determined by the ratio fc_/A:-|_. Away 
from any resonance, photon leaking favors the 0— photon 
(dark) state and thus, the corresponding interwell relax- 
ation is fast, i.e., k+ ^ fc_. Close to a resonance for 
A = 1/N (N integer), if the driving is resonant or the in- 
coherent resonant tunneling rate k- = flj^/Tiy is of the 
same order of fc+, the response is qualitatively modified 
and the resonant-antiresonant transition is now governed 
by the ratio fc_/fc+. 

Next, we complete our discussion with numerical re- 
sults for a concrete example. In Fig. [71 we plot (a) the 
overall oscillation amplitude, (b) the overall phase, (c) 
the rescaled quasienergy spectrum, and (d) the occupa- 
tion probabilities of the first three states j-i/'^) (n = 0, 1, 2) 
at the well bottom as a function of the effective Planck 
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FIG. 8. Solid line: Smallest eigenvalue of the Lindblad master 
equation. Dashed red line: the same without dissipative spin 
flips. Dashed-dotted orange line: decay rate for Jt/;*) — > \4'o)- 
Dotted lines correspond to k/2 and k/4. 



constant A. For A « 1/A^ {N = 2,3,4,5), we observe 
a resonant out-of-phase response (see Figs. [7] a) and b), 
whereas for A ~ 1/6, an antiresonance appears which 
is in phase with the modulation signal. For the same 
values of A, the occupation probabilities of the states in 
the well display peaks, see Fig. [7]d). The heights of the 
respective peaks are suppressed exponentially with N. 
According to the detailed balance condition in Eq. ([5T|). 
the ratio of the probabilities for nearby states close to the 
bottom of the well is exp(— /3cff). In fact, in logarithmic 
scale, the individual curves of the probabilities tends to 
be equidistant. The theoretical value for the gap between 
the curves is shown by the black double arrow, indicating 
that the agreement with Eq. pip is also quantitative. 

In order to underpin the drawn picture by more quan- 
titative results, it is is helpful to consider the eigenval- 
ues of the Liouville operator. For this, we start deep in 
the semiclassical regime, i.e., with small detuning. Here, 
a clear separation of time scales for the dissipative dy- 
namics on the bistable quasienergy surface occurs. Well- 
defined energy wells with a large quasienergy barrier in 
between exist for small detuning and allow for a clear 
description in terms of a single relaxation rate (see also 
Ref. |39| for a comprehensive review). In this regime, we 
find a single eigenvalue F which consists of the sum of fc_ 
and fc+ , is real and much smaller than the real parts of 
all the other eigenvalues. This is shown in Fig. [5] as black 
solid line. In order to emphasize the role of bath induced 
spin-flips, we have also computed the smallest eigenvalue 
F of the Liouvillian after removing those transition by 
hand from the master equation, see red dashed line in 
Fig. m For increasing A, i.e., when A > 0.25, we enter 
a regime, where the separation of time scales is not so 
clearly expressed and F becomes comparable to the real 
parts of the next three eigenvalues. The latter correspond 
to relaxation (one real eigenvalue) and to decoherence (a 
complex conjugated pair of eigenvalues), involving the 
pair of states \(j)i+) and \4>o)- We do not show them 
in the figure, but instead show the perturbatively deter- 



mined values k/2 (relaxation) and k/4 (decoherence) out 
of resonance as dotted horizontal lines. The peaks in the 
total rate F are due to resonant dynamical tunneling for 
fc+ as discussed in Section flVBI Due to the logarithmic 
scale, peaks in F are only visible when fc_)_ is of the order 
or larger than the background value given by fc_ . 



VII. CONNECTION TO OTHER MODELS 

The analysis presented here can in general be extended 
to any driven nonlinear oscillator coupled bilinearly to a 
thermal bath. For example, the Duffing oscillator is char- 
acterized by two classical stable solutions with opposite 
oscillation phase. In the quantum regime, two quantum 
squeezed states correspond to them. For weak driving, 
the small-oscillation solution can be identified with the 
0— photon quantum state. At low thermal energies, it is 
dominantly populated away from any resonance due to 
photons leaking into the bath. Thus, it can be regarded 
as stable in absence of any multiphoton transitions. How- 
ever, this stable quasienergy state is associated to a rela- 
tive quasienergy maximum, which is in direct contrast 
to the case of a static bistable potential. There, the 
lowest energy state is always the minimum of the true 
potential. At a multiphoton resonance, the zero-photon 
quasienergy state is no longer stable since excitations to 
a A^— photon state occur. Hence, it becomes metastable 
and generates a (anti-)resonance of the stationary oscilla- 
tion. This behavior of the quantum Duffing oscillator has 
been already predicted^i"— , but the link to the semiclassi- 
cal picture^i^ has not been drawn. Thus, the generic be- 
havior of a driven damped nonlinear quantum oscillator 
includes dynamically generated metastable states from 
which the system can escape via thermal diffusion, quan- 
tum activation or dynamical tunneling. In the regime of 
many photons in the resonator, the escape rate of dy- 
namical tunneling processes can be obtained in a semi- 
classical description, while in the regime of only few pho- 
tons (deep quantum regime), the escape can occur via 
resonant dynamical tunneling, leading to resonant mul- 
tiphoton transitions. Depending on the ratio of the Rabi 
frequency and the lifetime of the multiphoton state, the 
resonant dynamical tunneling can be coherent (for large 
quality factors) or incoherent (for small quality factors). 
Depending on the phase of the associated multiphoton 
transitions, a resonant or an antiresonant nonlinear re- 
sponse may arise. Such a situation is also expected in 
a Josephson bifurcation amplifier— operated in the deep 
quantum regime (note that all related experimental se- 
tups realized so far operate in the classical regime). 



VIII. CONCLUSIONS 

In conclusion, inspired by recent experiments, we have 
shown that a driven circuit QED setup can acquire a 
dynamical bistability. The relevant model to describe 
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this is the driven dissipativc Jaynes-Cummings model for 
which we have analyzed its nonlinear response properties. 
We have shown that a quasienergy surface can be de- 
rived in a rotating frame picture which clearly shows two 
metastable basins of attraction. This picture is also con- 
venient for studying the semiclassical limit. We have pre- 
dicted the existence of a metastable quantum squeezed 
state in the semiclassical limit and have discussed a con- 
nection between effective local temperature at the bot- 
tom of the quasienergy well and the squeezing parame- 
ter. We have analyzed the escape mechanisms from the 
metastable states and found resonant dynamical tunnel- 



ing, both in an incoherent and a coherent version. Our 
analysis adds another example to the series of nonlin- 
ear driven dissipative quantum systems with surprising 
and counterintuitive, but generic features. The thorough 
experimental investigation of these intrinsic quantum ef- 
fects is an interesting prospect. 
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